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Abstract 

In many SUSY models the first SUSY signal in the proposed International Linear 
Collider is expected to come from the pair production of fi, followed by its decay 
into r+ LSP. We study a simple and robust method of measuring the polarization 
of this r in its 1-prong hadronic decay channel, and show how it can be used to 
discriminate between SUSY models and to determine SUSY parameters. 



Introduction 



It is widely recognized now that the LHC should provide unambiguous signals of super- 
symmetry; but one needs the data from the proposed International Linear Collider(ILC) to 
discriminate between different SUSY models as well as for model independent determination 
of the SUSY parameters [1]. One expects the lighter tau slepton, 

fi — fRsmO-,- + Tl cos 9r, (1) 

to play a very important role in the SUSY study at ILC, since it is expected to be the 
lightest sfermion and the next to lightest superparticles (NLSP) over a large range of SUSY 
parameters in many important versions of the MSSM [2, 3]. Consequently the first SUSY 
signal at ILC is expected to come from the pair production of fi. 

In the universal SUGRA model the superparticle masses at the weak scale are related 
to the universal SUSY breaking mass parameters at the GUT scale, mo and mi/2, by the 
renormalization group equations(RGE). The lighter neutralinos(Zi 2) and chargino state(VFi) 
are dominated by the U(l) and SU(2) gauginos, B and W, with masses 



Ml = ~ O.5M2, 

00:2 

M2 = — mi/2 ~ 0.8mi/2. (2) 



The left and right slepton masses are given by 



mf^ = mo + 0.15mi/2- (3) 



At moderate to large tan/5(> 5) the tl^r masses are suppressed by negative contributions 
to the above RGE from the r yukawa coupling, 

which is twice as large in tr as in fi- Moreover there is significant mixing between the tl^r 
states, as represented by the off-diagonal term in the mass matrix 

= — m,-(^T- -|- /xtan/3) (5) 

which suppresses the mass of the lighter eigenstate f\ further down. Thus fi is predicted to 
be the lightest sfermion over the moderate to large tan/3(> 5) region, which is favored by 
the LEP limit on the light Higgs boson mass^. Moreover from eqs. 2-5 one expects fi to be 

-'^ Although there is no strict LEP limit on tan/3 for the maximal mixing scenario corresponding to At 
•\/6toj ~ 2.5 TeV, a moderate value of < 1 TeV implies small mixing and tan/3 > 5 [3, 4]. 
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the NLSP after Zi, over half the parameter space 



mo < mi/2. 



(6) 



The RGE of eqs.2 and 3 may not hold in non-universal versions of SUGRA and the gauge 
or anomaly mediated supersymmetry breaking models (GMSB and AMSB). Nonethless the 
suppressions coming from the Yukawa coupling and the off-diagonal mass term of eqs. 4 and 
5 continue to hold in these models for moderate to large values of tan /?. Consequently the 
fi is expected to be the NLSP over important ranges of parameter space in these models as 
well. For simplicity we shall consider the CP conserving versions of these models. 

The superparticle spectrum for a set of benchmark points for the above mentioned SUSY 
models, called the Snowmass points and slopes(SPS), are given in Ref. [5]. They satisfy 
all the experimental constraints including the cosmological constraint on the SUSY dark 
matter (DM) relic density. Three of the five universal SUGRA points have fi as the NLSP, 
including the typical SUGRA points(SPSl). The two exceptions represent the focus point 
region mo ^ 1111/2 and the so called funnel region(tan/3 = 55 and mo > mi/2). The 
DM constraint is satisfied through the pair annihilation of Zi by their coupling to the Z 
boson via their higgsino components in the focus point region and by their coupling to the 
pseudoscalar Higgs boson in the funnel region. In the other cases it is satisfied through the 
pair annihilation of Zi via a relatively light fi exchange ZiZi t^t^ or its co-annihilation 
with a nearly dengenerate fi. The list contains a non- universal SUGRA point (SPS6), where 
the Zi has a significant higgsino component and a roughly similar mass as Z2 and W^. In 
this case fi is slightly heavier than this cluster. It should be added here that specific non- 
universal SUGRA models have been considered, where the gaugino masses arise from the 
vacuum expectation value(vev) of a non-singlet chiral supcrficld, belonging to the 75 or 200 
representations of the GUT SU(5) group [6]. In these models the LSP consists of a set of 
degenerate higgsino Z-1^2 and W^. The fi turns out to be the next heavier state over the 
region(6) for moderate to large tan/3. The list of Ref.[5] contains two GMSB points, where 
the fi is the NLSP to a very light gravitino G in one case and it is slightly heavier than the 
Zi in the other. Finally it contains a AMSB point, where the fi is the NLSP to a degenerate 
pair of Zi, Wi LSP, representing the W . 

One sees from the above discussion that the fi is expected to be the NLSP in a wide class 
of SUSY models. Thus one expects an unambiguous SUSY signal from the pair production 
process 



In particular one can have a right polarized electron beam (Pg = -|-1), which couples only 
to the U{1)y gauge boson B. Thus for s ^ m|, (j{tr) = 4(T(fi); so that one can predict the 
fi pair production cross section in terms of its mass and mixing angle [7, 8]. The mass can 
be estimated from a threshold scan and the mixing angle from the production cross section. 




(7) 



followed by the decays 




(8) 
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Hence one can then predict the polarization of r coming from the decay (8) in terms of the 
composition of the LSP, i.e. 



W,f + (afi)" 

afi — ^Niitan9wsm9r 7=-^^ — Ni^cosOr, 

V2 \/2mw cos f3 

«n = 4k [^12 + ^11 tan 9w] cos 9^ - ^ ^""^ — -AT^gsin^^, (9) 
V2 V2mw cos p 



where 



Zi = NiiB + N,2W + NisHi + NuH2, (10) 

and we have made the coUinear approximation (mi- -C rrif^). 

In this work we investigate the prospect of using r polarization to probe the composition 
of Zi at the ILC. We shall use a simple and powerful method of measuring r polarization via 
its inclusive 1-prong hadronic decay, which was suggested in [9, 10] in the context of charged 
Higgs boson signal. The efficacy of this method for H"^ search has been corroborated by 
detailed simulation studies by both the CMS and ATLAS groups at LHC [11]. More recently 
the method has been used in the context of the fi signature at Tevatron and LHC [12]. 
However this seems to us to be the first application of this method for a e~^e~ collider signal. 
It should be noted here that some of the f signal studies at e+e~ collider have used the 
effect of T polarization in its exclusive decay channels r ^ nu [7, 8] and r ^ pu [7]. This 
has its own advantage if one wants e.g. to reconstruct the f mass from its decay kinematics 
rather than via threshold scan. But as a method of simply measuring the r polarization 
the inclusive channel has several advantages over the exclusive ones, since the latter suffers 
from lower statistics due to the branching fraction of the particular exclusive channel as well 
as larger systematic error due to the contaminations from the other decay channels. The 
estimation of this systematic error depends heavily on the ILC detector parameters, which 
is beyond the scope of this paper. We can only refer the interested reader to the paper of 
Nojiri et. al. [7] who have done a Monte Carlo simulation incorporating the parameters of a 
specific model detector. In particular one can see from their fig. 7 a significant contamination 
of the T pv decay channel from r aiv. On the other hand the method advocated 
here uses a feature of r polarization, to which the tt, p and ai decay channels contribute 
coherently. This obviates the necessity to separate these exclusive decay channels from one 
another as we sec below. In view of the expected improvements in a future ILC detector, it 
is quite possible that the systematic errors in the exclusive analysis will be reduced. In such 
a situation, the information obtained from both the inclusive and the exclusive analysis, can 
be used effectively together. 
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r Polarization 



The hadronic decay channel of r is known to be sensitive to r polarization [9, 10]. Wc shall 
concentrate on the 1-prong hadronic decay of r, which is best suited for r identification [11]. 
It accounts for 80% of r hadronic decay and 50% of its total decay width. The main 
contributors to the 1-prong hadronic decay are [3], 

T± ^ 7r±i/(12.5%), p±i/(26%), afu{7.5%), 

where the branching fractions for tt and p include the small K and K* contributions re- 
spectively, which have identical polarization effects. Together they account for 90% of the 
1-prong hadronic decay. The CM angular distribution of r decay into tt or a vector meson 
v{= p, oi) is simply related to its polarization via 

^ "^^^ ^:i+p.cos^) 



d COS 9 2 

T^dcose ~ ml + 2ml^ ^ ^ 

where L, T denote the longitudinal and transverse states of the vector meson. The fraction x 
of the r lab momentum carried by its decay meson is related in each case to 9 in the coUinear 
approximation via 

x=-{l + cos9) + - cos^). (12) 

The only measurable r momentum is the visible momentum of the r-jet, 

Pr-jet ^XPr. (13) 

It is clear from eqs. 11-13 that the hard r-jets are dominated by 71,Pl and an for Pr = +1, 
while they are dominated by pr and air for = —1. The two can be distinguished by 
exploiting the fact that the transverse p and oi decays favor even sharing of momentum 
among the decay pions, while the longitudinal p and Oi decays favor uneven sharing, where 
the charged pion carries either very little or most of the vector meson momentum. Therefore 
the fraction of the visible r-jet momentum carried by the charged prong, 

R = Pn^/Pr-jet: (14) 

is predicted to be peaked at the two ends i? < 0.2 and R > 0.8 for P^- = -|-1, while it is 
peaked at the middle for P^ = — 1 [9, 10]. Thus the R distribution of the hard r-jet can 
be used effectively to measure P^-, where the numerator can be measured from the tracker 
and the denominator from the calorimetric energy deposit of the r-jet. We shall use a jet 
hardness cut of 

> 25 GeV, cos 9r-jet < 0.75 (15) 
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for r-identification. But for the sake of comparison we shall show the results for a harder cut 
of p!^_jet > 50 GeV as well. It is well known that this signal can be effectively separated from 
the standard model background via modest cuts on the acoplanarity of the two r-jets and 
the missing-pT[7]. However we shall not consider these cuts here, since they do not affect 
the measurement of Pr- 

In our Monte Carlo simulation we shall use the TAUOLA packages [13] for polarized r decay 
into inclusive 1-prong hadronic channel. In addition to the above mentioned resonances 
it includes the small nonresonant contribution to the 1-prong hadronic r decay as well. 
There is some uncertainty in the paramctcrizations of the 7r,p, ai and the non-rcsonant 
contributions which have been tuned to the r decay data of different experiments by the 
respective collaborations. Indeed we expect this to be the main source of uncertainty in 
estimating the r-polarization from the R-distribution. To estimate the rough size of this 
uncertainty we shall compare the results obtained with the parameterizations of the ALEPH 
and the CLEO collaborations [13]. Wc shall use the PYTHIA event generator [14] to simulate 
the fi pair production and decay(eqs. 7, 8) at the ILC. 

For illustrative purpose we shall show the R distribution of hard r-jets coming from (7) 
and (8) for 

= 350 GeV, ruf^ = 150 GeV, m^^ = 100 GeV, (16) 

and different values of P-^, corresponding to different SUSY models. These masses are close 
to the fi and the Zi masses of the typical SUGRA point(SPSl) as well as the other relevant 
points of Ref. [5]. The size of the signal cross section for this beam energy and the fi mass is 
~ 200/6, while the typical luminosity for ILC is ~ 100/6~^. Correspondingly we have done 
a Monte Carlo simulation of 2x10^ signal events. 

We shall consider three different strategies for using P^ to probe the SUSY model and 
determine the model parameters: (A) Discrimination between SUSY models, (B) Partial de- 
termination of SUSY parameters and (C) Complete determination of EW SUSY parameters 
in association with the measurements of chargino iWi) cross-section and mass. 

A) Discrimination Between SUSY models: 

We consider four different SUSY models which are different versions of the MSSM - i) 
Universal Supergravity, ii) Higgsino LSP, iii) Anomaly Mediated SUSY Breaking(AMSB), 
iv) Gauge Mediate SUSY Breaking(GMSB). 

i) In the universal SUGRA model the LSP is dominated by the B component {Nn). 
Moreover at small tan/? the mixing angle cos^,- is small. Thus we see from eq.(9) that 

Pr ^ +1. (17) 

Since the mixing angle is determined by the off-diagonal term of eq. 5, large tan/? 
corresponds to large mixing (cos^^). But there is an effective cancellation between the 
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two terms of a^^ in eq.(9), so that Pt remains very close to +1 over practically the entire 
parameter space of interest [12]. Indeed it was shown in the second paper of Ref. [12] 
that Pr >0.9 throughout the allowed SUGRA parameter space, while Pr >0.95 in the 
mo < mi/2 region of our interest. Thus eq.l7 holds to a very good approximation at 
large tan (3 as well. 

ii) In nonuniversal SUGRA models, the gauge kinetic function and the resulting gaugino 
masses can be determined by a non-singlct chiral superficld, at the GUT scale, belong- 
ing to the SU{5) representations 24, 75 or 200. For the 75 and 200 representations the 
LSP is dominated by the Higgsino component over most of the parameter space [6]. 
Thus from eq. (9) 

~ cos^^^ -sin^^^, (18) 

particularly at large tan/?, where Higgsino couplings are enhanced. It should be men- 
tioned here that in this case there is a nearly degenerate pair of light neutrahnos 
Zi 2, which are both dominated by the higgsino components. Thus the predicted r 
polarizations of eq. 18 holds for each of them. 

iii) In AMSB model the LSP is dominated by the wino component (A'"i2) [15], i.e. 

Pr ~ -1. (19) 

iv) Finally in GMSB the LSP is the gravitino G, while the fi is expected to be the NLSP 
over a large part of parameter space [16]. In particular this is true for one of the two 
benchmark points of GMSB in Ref. [5] . Thus the fi — > tG decay implies 

P^ = sin2^^-cos^^^. (20) 

It may be added here that for the other GMSB benchmark point of [5] the fi is slightly 
heavier than a B dominated Zi, so that one expects essentially the same as in the 
universal SUGRA model. 

Thus for cos^i- = 1/2, which is a reasonable value of the mixing angle in the large tan /3 
region, we have 

Pr{i, a, iii, iv) ~ +1, -1/2, -1, +1/2. (21) 

Fig.l shows the resulting R distributions for the four P^ values corresponding to p^^j^t > 
25 GeV(upper pannel) andp!^_jgj > 50 GeV(lower panel). The last bin near R = 1 represents 
the — > n^u contribution. Thus the R distribution is expected to peak at the two ends for 
Pr >0 states while it is expected to peak near the middle along with a depleted — > n^u 
contribution for P^ <0. The difference increases with increasing Pr-jet cut as expected from 
eq. 11-14. But even for the p^_j^t >25 GcV cut of eq. 15, the differences in the R distribution 
are sufficient to distinguish the four models from one another. 
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Fig. 2 shows the fraction of events in the interval 0.2< R <0.8, 

^^.(0.2<ii<0.8) 

against the r polarization. This fraction is seen to decrease from 65% at = — 1 to 35% 
at Pr =+1 for the p^^j^t >25 GcV cut (solid lines); and the decrease becomes steeper 
for the harder p^^j^t > 50 GeV cut (dashed lines) as expected. Thus this fractional cross 
section can be used to measure Pr- In case the R < 0.2 region is inaccessible due to the 
difficulty in r identification for a soft charged track, the cross section a{0.2 < R) can be 
used for normalization in the denominator of eq. 22. To get a rough estimate of the size 
of uncertainty in the measurement of P^, wc have shown the predictions using the ALEPH 
and CLEO parameterizations of r decay [13] by two lines for each case. The horizontal 
spread between them corresponds to a APt- = ±0.03 (±0.05) near Pr=— 1(+1). One should 
of course add to this the uncertainty in the measurement of f (eq. 22). However being a 
ratio of measured cross-sections this quantity is dominated by the statistical error, which is 
~0.01 corresponding to ~ 10^ signal events. Hence the dominant error in P^ measurement 
is expected to come from the parametrization of r decay, as represented by the horizontal 
spread between the two lines of fig. 2. It should be noted here that even a conservative 
estimate of AP^ = ±0.1 would imply that the above four models can be distinguished from 
one another at the 5a level. 



B) Partial Determination of SUSY Parameters: 

In this section we shall assume the mass relation of eq.2 between the two electroweak gaug- 
inos, which holds in a fairly broad class of SUSY models. It holds in most SUSY models 
satisfying GUT symmetry, in which case the gaugino masses arise from the vev of a GUT 
singlet chiral superfield. This includes the SUGRA models with nonuniversal scalar masses. 
It holds for the GMSB as well. We shall assume that the tan /3 parameter can be indepen- 
dently measured from the Higgs sector at LHC or ILC. Moreover we shall assume that the 
absolute scale of one of the SUSY masses, representing some combination of Mi, M2 and 
II can also be measured independently. Then the relative magnitude of the Mi and the 
parameters can be estimated from the r Polarization. 

We have calculated the composition of Zi, as a function of the ratio Mi/|/x|, assuming 
the Zi mass of as in eq.l6 to fix the SUSY mass scale. Fig. 3 shows the resulting Pr from 
eqs.(9,10) for two representative values of tan (3 = 40 and 10 and cos^^f = 0.5 and 0.2. It 
shows that Pr is insensitive to this ratio for Mi/|/x| <0.5 at tan/3=40 (b,d) and Mi/|/i| <1 
at tan/3=10(a,c). But above these limits one can measure this ratio to an accuracy of ~ 
10% with a APr =0.1. 
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C) Complete determination of EW SUSY parameters 
in association with chargino (Wi) mass and cross section 
measurements: 

In most of the SUSY benchmark points of Ref. [5] discussed above, the Ughter chargino 
Wi mass lies in the range of 200-250 GeV. Therefore one hopes to see pair production of 
chargino in the first phase of ILC with a CM energy of ~500 GeV. In that case one can 
estimate some of the above mentioned SUSY parameters from this process, as discussed in 
[17]. In particular the mass and the composition of Wi can be easily obtained in terms of 
M2, IJ, and tan/3 by diagonalizing the chargino mass matrix. 

Thus one knows the couphng of a Wi pair to the Z boson in terms of these three parameters, 
while it has the universal EM coupling to 7. With a right polarized electron beam the Wi 
pair production proceeds only through s-channel 7 and Z exchanges. Thus by measuring 
this production cross section ctr along with the Wi mass from a threshold scan, one can 
estimate the composition of Wi. These two measurements determine two of the above three 
parameters. Hopefully the third parameter (tan (3) can be independently estimated from the 
Higgs sector at the LHC or the ILC, as mentioned above. Alternatively, one can measure 
the Wi pair production cross sections with left and transversally polarized electron beams, 
ai and (Tt- These two cross sections get contributions from t-channel z/g exchange along 
with the s-channel 7 and Z exchanges. Thus they depend on Ug mass along with the above 
three parameters. Nonetheless one can combine the three cross sections a 1,(^11 and cxr along 
with the Wi mass measurement to determine all the four parameters, as discussed in [17]. 
That would leave one SUSY parameter from the EW sector, (Mi), which cannot be directly 
measured at ILC in the absence of a Zi pair production signal. However one can use instead 
the Pr from the fi tZ\ decay to determine this parameter. Recall that along with Mi 
P-r has a dependence also on the stau sector mixing angle and tan/3. As a matter of fact, 
the role of Pr to probe the Higgsino component in the Zi, especially at large tan/3, has 
been emphasized [8]. As mentioned above, in the discussion we are assuming that these two 
parameters may be determined from other sectors. Of course, in a global fit, Pr will be one 
of the important variables in a model independent determination of Qr-, M2, Mi, /i and tan/3. 

Since both M2 and would be known from chargino production, the r polarization can 
be used to determine Mi relative to either one of them. For illustrative purpose we have 
assumed M2 = 250 GeV and three values of the ratio \\x\lM2 —2, 1 and 0.5. As in case 
of Fig.3 in this illustration also, we take two widely different sets of tan /3, cos ^^t-, chosen to 
span the range of expectations for these two, to demonstrate the dependence of Pr on them. 
The resulting polarization(PT-) is shown against M1/M2 in Fig. 4 for tan/3 = 40, cos ^^=0.5 
and in Fig. 5 for tan/3=10, cos ^,-=0.2. One sees from these figures that while Pr is fiat for 
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extreme values of M1/M2, it is a sensitive function of this ratio for 0.5 - M1/M2 ~ 2. One 
can measure the ratio M1/M2 to an accuracy of 5-10% over this interval with a AP^ =0.1. 
Thus knowing both n and M2 from Wi production one can use Pr to determine Mi relative 
to M2. In other words one can combine the Wi mass and cross section measurements with 
the r polarization to make a complete determination of the EW SUSY parameters. 

It should be mentioned here that there is an alternative method of determining Mi 
from kinematic distributions as suggested in ref. [18]. We are not in a position to make a 
comparative assessment of the relative merits of the two methods. Instead we would like 
to emphasize their complementarity. The kinematic distributions probe the LSP(Zi) mass, 
while the Pr measurement probes its composition. In this respect Pr plays a role analogous to 
the production cross section of chargino (Wi). Note that one can combine the measurements 
of Wi mass and composition via threshold scan and the production cross section ctr with 
those of Zi mass and composition via the kinematic distribution method and Pr- Thereby 
one can determine all the four parameters (/x, Mi, M2 and tan f3) with only the right polarized 
electron beam. On the other hand one can combine these four measurements with ai and 
ax- Then the resulting redundancy of information can be used as a consistency check of the 
CP conserving MSSM. Alternatively, it can be used to extend the analysis and probe for CP 
violation in MSSM [19]. 
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Figure 1: The normalized distributions in the fraction of T-jet momentum carried by the 
charged track for Pt= +1, +1/2, -1/2 and -1. The upper and lower panels correspond to 
p^^jet >25 GeV and p^_jet >50 GeV respectiyjly. 
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Figure 2: The fractional cross sections in the interval 0.2 < R < 0.8 (eq. 22) plotted against 
T polarization for p^_jgf >25 GeV (sohd lines) and p^_jgt >50 GeV (dashed hnes). The 
upper and lower lines for each case correspond to the parameterizations of r decay by the 
ALEPH and the CLEO collaborations [13] respectively. 
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Figure 3: Polarization of r shown against Mi/|^| for a fixed value of = 100 GeV 
and Mi/M2=0.5 and different sets of parameters: (a)tan/3 = 10, cos^,- =0.5; (b)tan/9 = 
A0,cos9f =0.5; (c)tan/5 = 10, cos^^ =0.2; (d)tan/5 = AQ^cosOf =0.2. 
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Figure 4: Polarization of r shown against M1/M2 for a fixed value of tan (3 = 40, cos9r =0.5, 
M2 =250 GeV and l/il/Ma =0.5,1,2. 
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